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Abstract. Lattice Hamiltonians with on-site interaction W have W — solutions, 
that is, many-body singlet eigenstates without double occupation. In particular, 
W = pairs give a clue to understand the pairing force in repulsive Hubbard models. 
These eigenstates are found in systems with high enough symmetry, like the square, 
hexagonal or triangular lattices. By a general theorem, we propose a systematic way 
to construct all the W — pairs of a given Hamiltonian. We also introduce a canonical 
transformation to calculate the effective interaction between the particles of such pairs. 
In geometries appropriate for the Cu02 planes of cuprate superconductors, armchair 
Carbon nanotubes or Cobalt Oxides planes, the dressed pair becomes a bound state 
in a physically relevant range of parameters. We also show that W = pairs quantize 
the magnetic flux like superconducting pairs do. The pairing mechanism breaks down 
in the presence of strong distortions. The W = pairs are also the building blocks for 
the antiferromagnetic ground state of the half-filled Hubbard model at weak coupling. 
Our analytical results for the 4x4 Hubbard square lattice, compared to available 
numerical data, demonstrate that the method, besides providing intuitive grasp on 
pairing, also has quantitative predictive power. We also consider including phonon 
effects in this scenario. Preliminary calculations with small clusters indicate that vector 
phonons hinder pairing while half-breathing modes are synergic with the W = pairing 
mechanism both at weak coupling and in the polaronic regime. 
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1. Introduction 



In the last decades, much effort has been devoted to exotic mechanisms of pairing and to 
the possibility of non- conventional superconductivity in new materials. In a list of the 
most important and frequently discussed examples, one could mention the cuprates[T], 
organic superconductors including fullerenes|2] and carbon nanotubes 0, ruthenatesjl], 
and Na-Co oxides jH]- Among the most interesting possibilities, correlation effects have 
been invoked, and the Hubbard model has been increasingly popular to achieve a 
simplified picture. 

The Hubbard model is believed to exhibit various interesting phenomena including 
antiferromagnetism, ferrimagnetism, ferromagnetism, metal-insulator transitions, etc. 
It is the simplest Hamiltonian covering both aspects of a strongly correlated electron 
system (like the Cu02 planes of cuprates), namely the competition between band-like 
behaviour and the tendency to atomic-like localisation driven by the screened Coulomb 
repulsion. Besides, several authors believe that it can exhibit a superconducting phase 
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in a certain parameters regime. Despite its simplicity, the Hubbard Hamiltonian 
cannot be exactly solved in more than Id and a large variety of approaches has been 
proposed to study the superconducting correlations in the ground state and at finite 
temperature. Bickers and co-workers|?j were among the first to propose the d^i-yi wave 
superconductivity in 2d. A recent survey of the superconducting properties of the single- 
and multi-band Hubbard model can be found in Ref.jSj. 

Another popular class of models are that obtained from the Hubbard Hamiltonian 
in the strong coupling limit. A large on-site Coulomb repulsive energy U between 
carriers of opposite spins tends to reduce the double occupancy. A straightforward 
perturbation theory in the parameters t of the kinetic term leads to Heisenberg-like 
Hamiltonians for fermions propagating through an antiferromagnetic background with 
an exchange interaction J between neighbouring spins. These are the so called t-J 
mo dels [Hj, [ID]- In the t-J model the double-occupancy is forbidden by the so called 
Gutzwiller projection [TT], and there is no on-site repulsion. However, superconductivity 
is a delicate phenomenon and the exclusion of double occupied sites costs kinetic energy, 
so the t-J model might fail to give an appropriate description of the Hubbard model 
when U is comparable to t. 

In this review article, we illustrate the 'W = 0" pairing mechanism, in which 
symmetry is capable of cutting down the on-site repulsion from the outset, without 
any need of the Gutzwiller projection. This effect works at any U/t and it relies on 
a configuration mixing which entails the presence of degenerate one-body states at the 
Fermi level. Moreover, we provide evidence that this configuration mixing, when applied 
to the full many-body problem, can produce pairing when physical parameter values are 
used. 

The pairing force in the Hubbard model is induced by repulsive interactions, 
and recalls the Kohn-Luttinger[T^ pairing in the Jellium. They pointed out that 
any three-dimensional Fermi liquid undergoes a superconducting instability by Cooper 
pairs of parallel spins and very large relative angular momentum /. A simplified view 
of the Kohn-Luttinger effect is given by considering one particle of the pair as an 
external charge. Then, the screening gives rise to a long-range oscillatory potential 
(Friedel oscillations) due to the singularity of the longitudinal dielectric function at 
2kF; here is the Fermi wavevector. The strict reasoning exploits the fact that the 
Legendre expansion coefficients of any regular direct interaction between particles of 
opposite momentum drops off exponentially in /. On the other hand, the second- 
order contribution to the scattering amplitude falls as l/l^ and at least for odd / leads 
to an attractive interaction. In the modern renormalization group language |13j. the 
second-order correction is obtained by summing up the marginal scattering amplitudes 
of the isotropic Fermi liquid coming from the so-called Forward channels, including, 
for antiparallel spins, a spin-flip diagram. This scenario does not work in the two- 
dimensional Fermi liquid, but going beyond the second-order perturbation theory the 
Kohn-Luttinger effect is recovered |14j. 

The present mechanism works with singlet pairs and differs from Kohn-Luttinger 
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one in other important ways. In contrast with the homogeneous electron gas, the lattice 
structure gives very tight pairs and the W = Q mechanism displays very clearly in tiny 
clusters as well. On top of that, it is worth noticing that in high-Tc superconductors 
the size of the Cooper pairs is expected to be of the order of few lattice constants 
and hence the pairing mechanism should lend itself to cluster studies. Macroscopically 
large lattices with periodic boundary conditions have large symmetry groups including 
the space groups (point symmetry + translational symmetry); in such conditions, the 
W = Q pairing mechanism is at its best. Similarly, in finite geometries, the largest 
binding energy is obtained in fully symmetric clusters while static distortions tend to 
unbind the pair. Exact calculations on finite models should bring to light interesting 
aspects of the microscopic origin of the pairing mechanism, and be useful as tests for 
the analytic developments. 

To test the superconducting nature of the pairs arising from repulsive interactions, 
one can use finite systems in gedanken experiments. We probe the behaviour of 
W = pairs in the presence of a static magnetic field and we show that they 
produce diamagnetic supercurrents that screen the vector potential. As a result the 
superconducting flux quantization is observed in various geometries. 

The paper is organised as follows. In Section we review the pairing mechanism in 
small symmetric clusters. We prove the existence of two-body singlet eigenstates with 
vanishing on-site Hubbard repulsion, that we call W = pairs. We have collected 
in Subsection 12.21 the somewhat more technical aspects, which are central for the 
mathematical foundation of the theory, while readers who are only interested in the 
phenomenology might skip it. We prove a general theorem on the allowed symmetries 
of such pairs. From the theorem we extract a practical recipe to build W = pairs in any 
symmetric geometry (finite or macroscopically large). A careful analysis on the smallest 
"allowed" cluster, the CUO4, shows that pairing can be obtained in a physical parameter 
range. The underlying pairing mechanism is investigated using many-body perturbation 
theory. In Section 01 we generalise the theory to arbitrary large systems. We introduce 
a non-perturbative canonical transformation leading to an effective Hamiltonian for the 
pair. The method is free from the limitations of perturbation theory; the relation of 
the present formalism to Cooper theory from one side and to cluster results from the 
other is discussed. Two kinds of bound states of different symmetries result, and the 
dependence of the binding energy on the filling and other parameters is explored. In 
Section 0] we study the Hubbard Model at half filling. We remove the ground-state 
degeneracy in first order perturbation theory by means of a suitable local formalism. 
We show that the ground state is the spin singlet projection of a determinantal state 
exhibiting the antiferromagnetic property: the translation by a lattice step is equivalent 
to a spin flip. As an illustration, the 4x4 square lattice is studied in detail. The half 
filled antiferromagnetic ground state is doped with two holes and an effective interaction 
between them is derived by means of the canonical transformation. The analytical 
results agree well with the numerical ones and this shows the predictive power of the 
approach. In SectionElwe investigate the 1^ = pairing mechanism in carbon nanotubes 
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and triangular cobalt oxides. Section IHl is aimed to study the superconducting magnetic 
response of symmetric Hubbard models with W = bound pairs. Section [7| deals with 
the inclusion of the lattice degrees of freedom; we show that phonons give a synergic 
contribution to the purely electronic mechanism and catastrophic Jahn- Teller distorsions 
do not occur. Finally, we draw our conclusions and outlook in Section |H1 

2. Pairing in the Hubbard Model 

In this review, we shall deal with Hubbard-like models of various geometries, designed for 
application to superconducting strongly correlated materials. The prototype Hubbard 
Hamiltonian reads 

H = K + W = tY. c],c,, + f/^n,tn,^ , (1) 

where K stands for the kinetic energy while W accounts for the on-site repulsive 
interaction. The summation on (z, j) runs over sites i and j which are nearest neighbours. 

In Eq. (P) the interaction term is repulsive and there is no electron-phonon coupling, 
so the very existence of pairing is a paradox. An effective attractive force comes 
from, e.g., the exchange of spin fluctuations but it must be stronger than the direct 
Hubbard repulsion to give rise a bound pair. Bickers and co-workers[Zj explored the 
consequences of a spin density wave instability on such a pairing force within the RPA 
approximation. They found a superconducting phase with pair- wavef unctions of dx2_y2 
symmetry in the 2d Hubbard model. The next level of calculations were carried out by 
using the FLEX approximation |il5j. This method treats the fluctuations in the magnetic, 
density and pairing channels in a self-consistent and conserving way. It was found that 
the antiferromagnetic fluctuations lead to a superconducting phase of dj.2_y2 symmetry 
which neighbours the SDW phase, in accordance with the previous findings. 

The phase diagram becomes less clear close to half filling because of the numerous 
infrared divergencies due to the nesting of the Fermi surface and Van Hove singularities. 
As a consequence, the results of any many-body treatment depend on the choice of 
diagrams to be summed. Renormalization Group (RG) methods^^l are a well controlled 
alternative approach to deal with Fermi systems having competing singularities. The 
RG has been used by several authors [TE],jT7j,[TH] to study the coupling flows at different 
particle densities. In agreement with the previous findings, RG calculations show a 
d-wave superconducting instability away from half filling. The underlying physical 
mechanism, namely exchange of spin- or charge-density fluctuations, is also the same as 
in the FLEX approach. 

Despite such evidence, there is no general agreement on the existence of a 
superconducting phase in the Hubbard modeljinj,[2n],j2I]- In order to clarify such a 
controversy it is very useful to have exact data to rely on. Unfortunately, the one 
band Hubbard model is exactly solved only in one spatial dimension [22] and no sign of 
superconductivity was found there jSHj- Very few exact analytical results are available 
for the 2 dimensional case [21]. Therefore in 2 dimensions the low- lying states must be 
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explored by means of exact diagonalizations on finite clusters. This kind of numerical 
calculation, complemented by analytical work and physical insight, may bring to light 
interesting local aspects of the microscopic pairing mechanism. 

The exact ground state of several small Hubbard systems have been numerically 
found by means of Lanczos and quantum Monte Carlo techniques. A pairing criterion 
that we shall discuss below was given by Richardson in the contest of nuclear physics [211 ■ 
Defining 

K{N + 2) = E{N + 2) + E{N) ~ 2E{N + I) (2) 

where E{N) is the A^-body ground state energy, 

A(A^ + 2) < (3) 

signals a bound pair in the ground state with + 2 particles and |A| is interpreted as 
the binding energy of the pair. Pairing was found under various conditions by several 
authors (SO], |SI], and among the geometries considered the 4x4 system [S2] is one of the 
most relevant for studying the pairing instability close to the antiferromagnetic phase. 

Below, we derive an analytic theory of pairing interactions, and report several 
detailed case studies where the formulas are validated by comparison with the numerical 
data. The results support general qualitative criteria for pairing induced by the on-site 
repulsion only. This shows that our approach successfully predicts the formation of 
bound pairs, and explains why other ingredients like strong off-site interactions are 
needed in other geometries. The analytic approach is also needed to predict what 
happens for large systems and in the thermodynamic limit. Increasing the cluster 
size the computed pair binding energies show a rapid decrease, and several authors 
on the basis of the numerical data consider pairing in the Hubbard model as a size 
effect. Unfortunately, the number of configurations grows in a prohibitive way with 
the cluster sizelSH] and numerical data currently available on 4x4 or even 6x6 clusters 
cannot provide reliable extrapolations to the bulk limit. In Section we show that we 
understand the trend analytically very well, that much larger cells (at least 30x30) are 
needed to estimate the asymptotic behaviour and that we have reasons to believe that 
pairing with a reduced but substantial binding energy persists in the full plane. 



2.1. W = Pairing in Cu-0 Clusters 

In this Section, we illustrate the concept of = pairs and the way they become bound 
states, by using examples with a geometry relevant for the cuprates. Our starting point 
is the three-band Hubbard Hamiltonian 

H = K + W + W^s-site (4) 

where 

K = t {p]^dia + h.c.) + tpp p]aPj'<T + SdY + ^pY1 
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Table 1. Character table of the Civ syrametry Group. Here 1 denotes the identity, 
C2 the 180 degrees rotation, c['^\ C^~^ the counterclockwise and clockwise 90 degrees 
rotation, a^, Uy the reflection with respect to the y = and x = axis and (7+, 
the reflection with respect to the x = y and x = ~y diagonals. The Civ symmetry 
Group has 4 one-dimensional irreducible representations (irreps) Ai, A2, i?i , B2 and 1 
two-dimensional irrep E. In the last column it is shown how each of them transforms 
under Civ 

and 

Here, pj (di) are fermionic operators that destroy holes at the Oxygen (Copper) ions 
labelled j (i) and n = p^p ( = Sd) is the number operator, (ij) refers to pairs of nearest 
neighbors i (Copper) and j (Oxygen) sites. The hopping terms correspond to the 
hybridization between nearest neighbors Cu and O atoms, and are roughly proportional 
to the overlap between localized Wannier orbitals. 

The parameters Ud and Up are positive constants that represent the repulsion 
between holes when they are at the same Copper id) and Oxygen (p) orbitals, 
respectively. Upd has a similar meaning, i.e., it corresponds to the Coulombic repulsion 
when two holes occupy two adiacent Cu-0 sites. The on-site energies Sp and Sd refer 
to the occupied orbitals of Oxygen and Copper. In the strong coupling hmit, and 
with one particle per unit cell, this model reduces to the spin Heisenberg Model with a 
superexchange antiferromagnetic coupling P^.pT]. 

From a band structure calculation and by best fitting the results of ah initio 
calculations^! one can roughly estimate the actual values of the parameters in the 
Hamiltonian of Eq. (jlj. Our preferred set is (in eV): Sp — Ed = 3.5, t = 1.3, tpp = —0.65, 
Ud = 5.3, Up = 6 and, most probably, Upd < 1.2. 

As Cini and BalzarottiJHl] pointed out, highly symmetric clusters possess 2- holes 
singlet eigenstates of H which do not feel the on-site repulsion W; such eigenstates 
were called W = pairs and play a crucial role for pairing. In order to have W = 
solutions, the clusters must possess the full 6*41, (square) symmetry, and must be centered 
around a Cu site. Group arguments are central to our approach and we show the 
Characters of the 6*4^, Group in Table ^ The symmetry requirements are so stringent 
that clusters with such properties had not been studied previously. In particular, the 
present discussion does not apply to the geometries, like those examined by Hirsch et 
a/.jHS],|SSl and Balseiro et a/.jnO], which are forbidden from our viewpoint. The CU4O4 
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Table 2. One-body levels of the CUO4 cluster in units of t as a function of the 
dimensionless parameter r = tpp/t. 

geometry considered by Ogata and Shibaj^T] has the C^y symmetry, but lacks the central 
Cu, and therefore it is forbidden. Those built from degenerate one-body levels are not 
W = pairs; they feel the on-site repulsion Up on the Oxygen sites, and yield A < 
only if Up is very small [SI], |SH1- 

Note that we are using the same t for all Cu-0 bonds, while some authors use 
other conventions. For instance, F. C. Zhang and T. M. Ricejl2] use an alternating 
sign prescription for t, which may be obtained by changing the sign of all the O orbitals 
in the horizontal lines containing Cu ions. The two pictures are related by a gauge 
transformation, under which the orbital symmetry labels Ai and Bi are interchanged. 
Some of the symmetry related information is gauge dependent and unobservable, while 
some is physical (e.g. degeneracies are). 

The reason why symmetry is so basic for our mechanism is that the W = 
two-body eigenstates arise when two holes occupy degenerate one-body levels. When 
such degenerate states are partially filled in the many-body ground state, bound pairs 
can form. The interacting ground state must be described in terms of a many-body 
configuration mixing, but the presence of the W = pair imparts to the many-body 
state a special character; it can be described in terms of a bound pair moving on a 
closed-shell background. This is the main point that we want to make, and it is deeply 
related to the symmetry quantum numbers. 

_•' I ■. 

I 
I 

I 
I 
I 

Figure 1. The topology of the CUO4 cluster. 

The smallest square Cu-0 cluster, CUO4, see Fig. d is also the simplest where 
W = pairing occurs. The Hamiltonian reads: 

HcnO, = tJ2idiP^<.+PlM+tpp E P.W + f/(^l'^^f^ + E^Sf^lf)-(5) 

icr <ij>,cr i 

The one-body levels and their symmetry labels are reported in Table El Labelling the 
Oxygen atomic sites by the numbers 1, 2, 3, 4, the diagonalizing creation operator are 



O Oxygen 
O Copper 
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t _ ^ M t 

4ia = \ {pL - ^2<x + pL - pla 



2 

CAia(l) = ^2_^^ ("+4 + I'L + ]'2a + I'Sa + P\a 
4l-(2) = 4 + I'L + I'L + P\a + l^l) (6) 

where a and /3 depend on r = tpp/t as follows 



4 (±l±r2 + rV4TT^ 



±5 r ± 2 r3 + V4+72 + 2 r2 V4+72 



2 particles 3 particles 4 particles 

Figure 2. Many body states for the Cu04 cluster. For each number of particle, it has 
been reported the component of highest weigth of the ground state. 

Let us build a 4-hole state in the CUO4 cluster in the non- interacting limit, 
according to the aufbau principle, see Fig. 2. The first two holes go into a bonding 
level of Ai symmetry; this is a totally symmetric (^^i) pair. For negative r, the other 
two holes go into a non-bonding level of E{x,y) symmetry, which contains 4 spin- 

/ 4 \ 

orbital states. The Pauli principle allows 12/"^ different pair-states. The irrep 

multiplication table allows for labelling them according with their space symmetry: 
E ^ E = B2. It is also straightforward to verify that A2 is a spin- 

triplet, ^A2, while the remaining irreps are spin-singlets, ^Ai, ^Bi and ^i?2- From Eq. 
(jni) one readily realize that the B2 singlet operator 

is a Ty = pair (no double occupation). Note that ^-82 is the symmetry label of the pair 
wave function in the gauge we are using, and must not be confused with the symmetry 
of the order parameter. 

To first order in perturbation theory, the 4-body singlet state of B2 symmetry is 
degenerate with the A2 triplet; Hund's rule would have predicted a ground state. 
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However, the true ground state turns out to be singlet, for reasons that we shall study 
below. The numerical results on the CUO4 cluster show that A (4) is negative for 
> r > —0.5 and that its minimum value occurs at r = 0, when the non-bonding 
orbitals Bi and E become degenerate: the paired interacting ground state is also 
degenerate in this particular case. A symmetry analysis of this accidental degeneracy 
is postponed to the next Section. In Fig. 01 we plot A (4) for r = 0, — = 0, 
Upd = and Up = Ud = U. A (4) has a minimum at U ~ 5 t and it is negative when 




Figure 3. A(4) (in t units) as a function of U/t. The maximum binding occurs at 
U ^ U where A(4) « -0.042 t. For U > 34.77 t (not shown), A(4) becomes positive 
and pairing disappears. 

< U < 34.77 t. We emphasize that A (4) becomes positive for large values of U/t and 
hence pairing disappears in the strong coupling regime. In the present problem U must 
exceed several tens of times t before the asymptotic strong coupling regime sets in. A 
perturbation theory will strictly apply at weak coupling where the second derivative of 
the curve is negative. However, qualitatively a weak coupling approach is rewarding in 
all the physically interesting range of parameters. The sign of A depends on U and r 
and its magnitude is unlike any of the input parameters; below we show that this new 
energy scale comes out from an interference between electron-hole exchanges of different 
symmetries. 

On the other hand, at positive r's the Bi non-bonding level is pushed below the 
degenerate one and A(4) becomes large and positive (at r = +0.65, A(4) = 0.53 eV). 

Next, we discuss the dependence of A (4) on the other parameters ||34j. If we 
decrease Ep, A (4) decreases because this makes the system more polarizable. The Sp 
dependence when all the other parameter are kept fixed and Upd = is almost linear 
down to Ep = 0. According to Ref.^^ positive Upd values do not spoil the mechanism, 
and tend to be synergic with it. Indeed, values of Upd > 0.6 eV give negative A(4) 
values even for r = —0.65 eV (in the range 0.2-1.2 eV considered in Ref . |34j A(4) is a 
monotonically decreasing function of Upd) ■ Finally, we have numerically studied how the 
distortions effect A (4). We have found that any lowering of the symmetry is reflected 
by a corresponding increase of A(4)|34j. 

In the CUO4 cluster we need four holes to have a paired ground state. Hence, 
the total hole concentration is ph = 0.8. This value is too large by a factor of 2 with 
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respect to the experimentally observed p/^ ~ 0.4 of the optimally doped systems. It is 
important to realize that these undesirable features are peculiar of the prototype CUO4 
cluster, and already disappear in CU5O4, the next larger cluster of the same symmetry. 
In fact, four holes are still sufficient to reach degenerate states, but ph ~ 0.44 is much 
closer to the experimental value. We have performed numerical explorations in other 
fully symmetric clusters like CU5O4 and CU5O16 and we have found negative values of 
A (4), of the order of few meV, using physical parameters. In all the allowed clusters 
up to 21 atoms, the lowest one-hole level belongs to Ai symmetry, and the next E level 
yields the = pair. The interactions produce a non-degenerate ^-82 4-hole ground 
state having the same symmetry as the W = Q pair. The interested reader may see 
Refs.|lH].|H] for the details. Below, in Section El we show [13] that the A(4) < arises 
from an effective attractive interaction between the holes of the W = Q pair; the same 
interaction is repulsive for triplet pairs. 

2.2. Symmetry of the W = Q pairs: a general theorem 

Since the mechanism depends on symmetry in such a fundamental way, we must refine 
the Group Theory analysis. We discovered a powerful and elegant criterion jlH] to 
construct all the W = Q pair eigenstates on a given lattice A by using projection 
operators. 

Let Qq be the the symmetry Group of the non-interacting Hubbard Hamiltonian 
K = '}2{ij),aUj(^j(,Ci<T- We assume that no degeneracy between one-body eigenstates is 
accidental, hence Qq must contain enough operations to justify all such degeneracies. 
Let us label each one-body eigenstate oi H = K + W, W = J^ieAUiUi^nii, with an 
irreducible representation (irrep) of Qq. 

Definition. An irrep 77 is represented in the one-body spectrum of H if at least one 
of the one-body levels belongs to 77. 

Let S be the set of the irreps of Qo which are represented in the one-body spectrum 
of H. Let be a two-body eigenstate of the non-interacting Hamiltonian with spin 
= and P*-^-* the projection operator on the irrep rj. We wish to prove the 

W = Theorem. - Any nonvanishing projection of on an irrep not contained in 
£, is an eigenstate of H with no double occupancy. The singlet component of this state 
is aW = pair. Conversely, any pair belonging to an irrep represented in the spectrum 
must have non- vanishing W expectation value, see Fig. (j3)): 

r] ^ S ^WP^'^'>\'^) = 0. (8) 

Remark: Suppose we perform a gauge change in H such that Qq is preserved; clearly, 
a. W = pair goes to another W = pair. Thus, the theorem implies a distinction 
between symmetry types which is gauge-independent. 
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Proof: Let us consider a two-body state of opposite spins belonging to the irrep r/ 
of Qq: 

Then we have 

We define P*-''^ as the projection operator on the irrep r]. Since 

ieA ieA 

if P^'^^\i T,i i) = Vi G A, then = Vi G A. It is worth to note that this 

condition is satisfied if and only if P^'^^\i(y) = Vi G A, where \ia) = cl^\0). Now it is 
always possible to write \ia) in the form \ia) = J2r)e£ '^^^K^)\v^) where \t]„) is the one- 
body eigenstate of H with spin a belonging to the irrep rj. Hence, if r^' ^ £, P^'^'^ia) = 
and so P'-'^'^li T, U) = 0. Therefore, if |?/^('')) is a two-hole singlet eigenstate of the kinetic 
term and f] ^ S, then it is also an eigenstate of W with vanishing eigenvalue, that means 
a IV = pair. 
Q.E.D. 

The content of the theorem is schematically represented by Fig. (j3)). 




Figure 4. Schematic of the W ~ theorem. 



We already know that CUO4 is a good example of this theorem. Indeed, the irrep 
B2 of the ly = pair is not represented in the spectrum (see Table |2)); A2 is also not 
represented, but it yields no two-body states. This possibility is admitted by Eq. (jS)). 

The particular case of CUO4 with r = is of special interest as an illustration. There 
is an accidental degeneracy between E{x, y) and Bi orbitals and therefore a three-times 
degenerate one-body level exists. With 4 interacting fermions, pairing occurs in two 
ways, as Ai and B2 are both degenerate ground states. This agrees with the theorem; 
in fact, the extra degeneracy cannot be explained in terms of C4^y, whose irreps have 
dimension 2 at most. For r = any permutation of the four Oxygen sites is actually a 
symmetry and therefore Qq is enlarged to 5*4 (the group of permutations of four objects). 
5*4 has the irreducible representations Ai (total-symmetric), B2 (total-antisymmetric), 
S (self-dual), 7i and its dual T2, of dimensions 1, 1, 2, 3 and 3, respectively. These 
irreps break in 6*41, as follows 



Ai = Ai, Ti = Bi® E, T2 = A2® E, B2 = B2, S = Ai® B2 
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Labelling the one-body levels with the irreps of 5*4, one finds that 8 is not contained in 
the spectrum and thus yields W = Q pairs: 

= Ml © . (9) 

The explicit expression of the corresponding pair-creation operators is 

,-|- 2-|--|- 

for the total-symmetric pair and Eq. ((7j) for the B2 component. 

In Section 14.31 we shall show use the theorem with a Qq that includes symmetries 
which are not isometries. We shall see that this theorem puts useful restrictions on the 
many-body ground state symmetry. The complete characterization of the symmetry of 
= pairs requires the knowledge of ^o- A partial use of the theorem is possible if 
one does not know Qq but knows a Subgroup Qq. It is then still granted that any pair 
belonging to an irrep of Qq not represented in the spectrum has the ly = property. 
On the other hand, accidental degeneracies occur with a Subgroup of Qq, and by mixing 
degenerate pairs belonging to irreps represented in the spectrum one can find W = Q 
pairs also there. 

The theorem tells us that the bigger is Qq, the larger is the number of = 
pairs. Indeed, for a given system, the number of one-body eigenvectors is fixed, while 
the number and the dimension of the irreps grow with the order of the group. Therefore, 
also the number of irreps not represented in £ grows, meaning more W = Q pairs. At the 
same time, a big "non-interacting" symmetry Group Qo grants large level degeneracies. 
As we have seen, this fact allows the existence oi W = Q pairs formed by degenerate 
orbitals. In the next Section we show that an anomalously low effective repulsion takes 
place in the interacting system when such levels are on the Fermi surface and that in a 
certain region of the parameters space this leads to pairing. 

3. Theory of Pairing in repulsive Hubbard Models 

In Section |21 we have shown that W = Q pairs may lead to a paired ground state in 
fully symmetric clusters centered around a Cu ion. In this Section we extend the 
theory to the full plane, and show that such pairs provide the natural explanation of the 
pairing instability of the Hubbard Model. By a novel canonical tranformation approach, 
we shall calculate the effective interaction W^q between two holes added to the ground 
state of the repulsive Hubbard Model. Furthermore, we shall show that W^s is attractive 
in the IV = pair-symmetry channels. The method is a particularly efficient way to 
perform the configuration interaction calculation. It is based on the symmetry and sheds 
light on the origin of pairing in the CUO2 plane. It is worth noticing that W = Q pairs 
emerge from symmetry alone and hence remain W = Q for any coupling strength U . 
This is the reason why weak coupling expansions provide often good approximations at 
intermediate coupling as well, as observed by several authors|TF] . jiH] . [1^ . 



CONTENTS 



14 



3.1. W = Q Pairs in the Planar Hubbard Model 

The Cu02 planes of cuprates have a large symmetry group which contains the space 
group Q = T^C^jj, where T is the group of translations and ® stands for the semidirect 
product. However, the Optimal Group Qq is much larger than ^jSO]- For this reason, 
instead of using the W = theorem, it is easier to follow the simple route of projecting a 
single determinantal state on the irreps of C^..,|51j.|3^. This partial use of the symmetry 
still gives enough solutions to demonstrate pairing. 

Let us focus on W = pairs with vanishing quasimomentum. Omitting band 
indices, we denote by 



a two-hole determinantal state, where k and — k label degenerate one-body eigenfun- 
ctions of K. The combination |(i[k]) + k]) is singlet and |(i[k]) — k]) is triplet. 
We introduce the determinants |(i[_Rk]) = |(i[kij]), R G C^^, and the projected states 



where x^^KR) is the character of the operation R in the irrep r]. 

In the Three Band Model we showed that W = pairs are obtained by projecting 
onto the irreps A2 and 521^3 (which are not represented in the one-body spectrum, as 
it should). These irreps differ by those obtained in the fully symmetric clusters, where 
W = pairs had Ai and B2 symmetry [23]. The reason for this change is a twofold 
size effect. On one hand, Ai pairs have the W = property only in clusters having the 
topology of a cross (whose symmetry group is S4 rather than C^^), but do not generalize 
as such to the full plane (when the symmetry is lowered to C^^). On the other hand, 
the small clusters admit no = solutions of A2 symmetry if only degenerate states 
are used. 

We recall that a necessary condition for pairing in clusters is that the least bound 
holes form a.W = pair, and this dictates conditions on the occupation number. In the 
full plane, however, W = pairs exist at the Fermi level for any filling. We also observe 
that the W = pairs obtained with the above procedure are not all the possible W = 
pairs in the One and Three Bands Hubbard models, since the W = theorem has not 
been fully exploited. However, it can be shown that they are the only W = pairs with 
holes in degenerate one-body levels and with vanishing quasimomentum. 

3.2. Canonical Transformation 

In this Section we intend to study the effective interaction Wes between the particles 
forming a. W = pair added to the A^-body interacting ground state \^u{N)). Since 
the two extra particles cannot interact directly (by definition of = pair), their 
effective interaction comes out from virtual electron-hole exchanges, and in principle 
can be attractive. 



\d[k]) = cl^cl^jO) 



(10) 




(11) 
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Many configurations contribute to the interacting (A^ + 2)-body ground state 
|$^(A^ + 2)) and we need a complete set S to expand it exactly; as long as it is complete, 
however, we can design S as we please. We can take the non-interacting A^-body Fermi 
sphere \^q{N)) as our vacuum and build the complete set in terms of excitations over 
it. In the subspace with vanishing spin z component, the simplest states that enter the 
configuration mixing are those obtained from |$o(-^)) by creating two extra holes |55j 
over it; we denote with \m) these states. Similarly, along with the pair m states, we 
introduce the 4-body a states, obtained from |<I>o(-^)) by creating 2 holes and 1 electron- 
hole (e-h) pair. Then, S includes the 6-body (3 states having 2 holes and 2 e-h pairs, and 
so on. We are using Greek indices for the configurations containing the electron-hole 
pairs, which here are playing largely the same role as phonons in the Cooper theory. By 
means of the complet set S we now expand the interacting ground state 

\^u{N + 2)) = ^a^|m) + ^a«|a) + Y.(^p\(3) + .... (12) 

m a 1^ 

and set up the Schrodinger equation 

H\^u{N + 2)) = E\^u{N + 2)). 

We emphasize that Eq. (|T^ is a configuration interaction, not a perturbative expansion. 
When the number N of holes in the system is such that |$o(-^)) is a single non- 
degenerate determinant (the Fermi surface is totally filled), the expansion (fT^ for the 
interacting ground state is unique and we can unambiguously define and calculate the 
effective interaction Wes- This is done by a canonical transformation . |51j . [52] of the 
many-body Hubbard Hamiltonian. We consider the effects of the operators K and W 
on the terms of |$t/(A^ + 2)). Choosing the m, a, P, . . . states to be eigenstates of the 
kinetic operator K, we have 

K\m) = E^\m), K\a) = E^la), K\(3) = Ep\(3) , .... 

The Hubbard interaction W can create or destroy up to 2 e-h pairs. Therefore, its action 
on an m state yields 

W\m) = YWm',n.\m') + + E ^(,,^(3), 

m' a /3 

on an a states yields 

W\a) = Y.W^,a\m) + E Wa',a\a') + + E ^7«l7), 

ma' P 1 

and so on. In this way we obtain an algebraic system for the coefficients of the 
configuration interaction of Eq. (fT^ . In order to test the instability of the Fermi 
liquid towards pairing, it is sufficient to study the amplitudes of the m states. In the 
weak coupling limit this can be done by truncating the expansion (fT^ to the a states. 
As we have shown the inclusion of the 7,... states produces a i?-dependent 
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renormalization of the matrix elements of higher order in VT, leaving the structure of 
the equations unaltered. 

By taking a linear combination of the a states in such a way that 

{K + = 5^^,E'^ 

the algebraic system reduces to 

[Em - E]am+ Yy^rn,m'am' + Xl^"^.""" = ^ (^3) 

m' Of 

[E; - a, + Y.^o.,m'am' = 0. 

m' 

Solving for and substituting back in Eq. (fT^ . we end up with an equation for the 
dressed pair lip) = J2m.'^m\^)- The effective Schrodinger equation for the pair reads 

{K + W + S[E]) 1^) = ifpairlV) = E\tlj) (14) 

where 

(5[^]k^' = -E%%^- (15) 

is the scattering operator. The matrix elements Wm,m' in Eq- (fT^ niay be written as 
the sum of two terms representing the direct interaction W^^^, between the particles 
forming the pair and the first-order self-energy Wm'- 

Similarly, in S[E] we may recognize two different contributions; one is the true effective 
interaction Wes, while the other is the forward scattering term E 

The first-order self-energy and the forward scattering term are diagonal in the indices 
m and m' . Wm and E^ renormalize the non-interacting energy E^ of the m states: 

Em ^ E(f) = Em + Wm + Em. 

Eq. fll4|) has the form of a Schrodinger equation with eigenvalue E for the added 
pair with the interaction M/^*^) -|- W^e- Here, the W = Q pairs are special because 
W^'^'^ vanishes. We interpret as the wave function of the dressed pair, which is 
acted upon by an effective Hamiltonian impair- This way of looking at Eq. (fT^ is 
perfectly consistent, despite the presence of the many-body eigenvalue E. Indeed, if 
the interaction is attractive and produces bound states, the spectrum contains discrete 
states below the threshold of the continuum (two-electron Fermi energy). This is a 
clear-cut criterion for pairing, which is exact in principle. The threshold is given by 

4^)^min[i?Cf)(E)], 

{m] 
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and contains all the pairwise interactions except those between the particles in the 
pair, i?^^ must be determined once Eq. p4|) has been solved (since F depends on 
the solution). The ground state energy E may be conveniently written as Ep + A. 
A < indicates a Cooper-like instability of the normal Fermi liquid and its magnitude 
represents the binding energy of the pair. 

It is worth noticing that in principle the canonical transformation is exact and it is 
not limited to the weak coupling regime. In the numerical calculations, however, some 
approximation is needed. In practice, we shall compute the bare quantities, that is, we 
shall neglect the 6-body and higher excitations in the calculation of and F. This is a 
resonable approximation if we compute small corrections to a Fermi liquid background, 
and the exact numerical results obtained in small clusters suggest that this is the case, 
see next Section. 

3.3. Pairing in the CuO^ Cluster 

We defined the pairing energy for a cluster with + 2 particles by introducing the 
quantity K{N + 2) = E{N + 2)+E{N)-2E{N + 1). At weak coupling we have shown in 
Ref . that A agrees with the binding energy A obtained from many-body perturbation 
theory. Here, we consider again the CUO4 cluster; we intend to test our canonical 
transformation by comparing the analytic results with the exact diagonalization data 
obtained in Section |21 

The CUO4 cluster can host a, W = Q pair of B2 symmetry. At weak coupling, 
we can ignore the effect of the renormalizations in Eq. (fT3j) . Furthermore, we can 
discard the inter-shell interactions. Indeed, the one-body levels of finite clusters are 
widely separated and the dominant contribution to the effective interation comes from 
the intra-shell interaction. Accordingly, in Eq. ()15|) we set the m = m' = B2. The 
one-body levels and their symmetry labels are reported in Table El 

The effective interaction in the B2 symmetry channel is 

2 1 

£bi + £Ai — E eA[ + £Ai — E _ 

The binding energy is obtained by expressing the lowest eigenvalue E of impair as 
E = Ep + A. The renormalized two-hole Fermi energy Ep coincides with the bare 
one at this level of accuracy, i.e., Ef^ = 2eA^. The results are shown in Fig. El for r = 0. 
We can see that A is negative and hence the renormalization approach predicts pairing, 
at least at weak coupling. We also observe that the order of magnitude of |A| is 10~H, 
which is not comparable to any of the U and t input parameters. The reason is that the 
interaction, which vanishes identically for the bare W = pairs, remains dynamically 
small for the dressed quasiparticles. This suggests that a weak coupling theory may be 
useful to study the pairing force, despite the fact that U is not small compared to t. 

Next, we intend to compare A with the quantity A (4) obtained from exact 
diagonalizations, see Fig. O At weak coupling the agreement is excellent. However, 
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Figure 5. Trend of A versus U/t in units of t. 

|A| is ~ 2 times greater than |A(4)| for U/t ^ 1. This means that the inter-shell 
interactions and the renormalizations of the a-state energies have an important weight 
in determining the right value of A. However, what is comfortable is that the analytic 
approach predicts the right trend of the binding energy. In the next Section we shall 
use the canonical transformation to study larger and more physical systems. 

3.4- Pairing in the Cu02 Plane: Numerieal Results 

Using the analytic expression for the effective interaction in the full plane, we have 
performed exploratory numerical estimates of A by working on supercells of A^a x 
cells, with periodic boundary conditions. For the sake of simplicity, we have neglected 
the minor contributions from the higher bands and considered the dominant intra-band 
processes, in which empty states belong to the bonding band. We have solved the 
problem in a virtually exact way for N\ up to 40. Several supercell calculations had 
been reported to date^Qj, but no conclusive evidence of a pairing instability was reached 
due to the difficulty of dealing with size effects. 

First, we have investigated triplet pairing but, as in the clusters, we have found a 
repulsive effective interaction. On the contrary, W = Q singlets show pairing, in line 
with our previous findings in small clusters [31], 5 Since screening excitations are 
explicitly accounted for in the Hamiltonian, it is likely that [/ is a bare (unscreened) 
quantity, which justify large values. A stronger interaction causes smaller pairs and 
speeds up convergence within attainable supercell sizes. In Table El we report the 
results for ^i?2 pairs with Sp = —1.3 eV [half filling corresponds to ep = —1.384 eV and 
we have used as input data the set of current parameters already used for clusters, that 
is (in eV) t = 1.3, Ep = 3.5, Ed = 0, Up = 12.5, Ud = 11.2]. 

With supercell sizes TVa > 40 calculations become hard. Thus, we need a simple 
solvable model in supercells and in the infinite plane to make reliable extrapolations 
of numerical results. To this end, we introduce an Average Effective Interaction (AEI) 
—Vcs, Vcs > 0, which is constant for all the empty states k and k'. For any A^a, is 
implicitly defined from 

_8_ 1 ^ ejE^-Ep) 

Kff Nl\'2{E^-Ep) + \Ay 

(the factor 8 comes from the projection onto the irrep B2). Although |A| decreases 
monotonically by increasing A^a, Kff remains fairly stable around 6 7 eV, see Table 01 
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p = NlNl 


-A (meV) 


Kff (eV) 




18 


1.13 


121.9 


7.8 


41.6 


20 


1.16 


42.2 


5.0 


9.0 


24 


1.14 


59.7 


7.0 


28.9 


30 


1.14 


56.0 


5.7 


13.2 


40 


1.16 


30.5 


6.6 


23.4 



Table 3. Binding energy of ^i?2 pairs in supercells. 



The relatively mild A^'a dependence of V^g supports the use of the AEI to extrapolate 
the results to the thermodynamic limit. The asymptotic value limTVA^oo |^| = Aasympt 
can then be obtained from 

8 _ /-o dep{e) 

Kff -Jep 2 (e — Sp) + l^asymptl ' 

where p is the density of states. 

The results for the pairs are seen to lead to bound states as well, with 
comparable A valueslHl]; the trend with doping is opposite, however, and the binding 
energy is nearly closing at Sp = —1.1 eV. 

Although the three-band Hubbard Model is an idealization of the strongly 
correlated Cu02 planes, it is interesting to observe that evidence of mixed {s + id) 
symmetry for the pairing state has been amply reported in angle-resolved photoemission 
studies [FT]. 

4. The Doped Hubbard Antiferromagnet 

The canonical transformation described in Section 13.21 relies on the uniqueness of the 
non- interacting ground state |$o(^))- l'^o(^)) is certainly unique if the Fermi surface 
is totally filled and it can be written as a Slater determinant. However, we are also 
particularly interested in the doped Hubbard antiferromagnet, and the antiferromagnetic 
ground state occurs at half filling, not in a filled-shell situation. 

We want to study the doped antiferromagnet since there are strong indications 
that the Fermi liquid is unstable towards pairing near half filling; they come from 
diagrammatic approaches |15j. renormalization group techniques|16j.jl7j and also cluster 
diagonalizations|321- Therefore exact results on the half filled Hubbard Model may be 
relevant to antiferromagnet ism and to the mechanism of the superconducting instability 
as well. The Lieb theorempS^ on the ground state spin-degeneracy of the half-filled 
Hubbard model ensures that for any bipartite lattice = AVJ with |^| = \B\, the 
ground state is unique for any interaction strenght U . Thus, we can use the canonical 
transformation with |$o(-^)) = lini(7^o+ \^u{N)) provided the number of holes A^ equals 
the number of sites |A|. Remarkably, the W = Q pairs are also essential tools to deal 
with the antiferromagnetic ground state solution. Below, we present two exact results. 
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First, we obtain the ground state |$^(|A|)) on the so called complete bipartite graph 
(CBG) for arbitrary but finite U . Then, we consider the square lattice and we resolve 
the degeneracy of the non-interacting ground-state multiplet at half filling. This latter 
result will then be used to study the pairing mechanism in the 4x4 Hubbard model. 

4-1. Half-Filled Hubbard Model on the Complete Bipartite Graph 

The Complete Bipartite Graph (CBG) A = ^ U ;B has bonds connecting any element 
of A with all the elements of B. In Fig. El we have drawn a few examples of finite- 
size systems. For = 1 and N = 2 the model is equivalent to a one dimensional 
ring of length L = 2, 4 respectively. For A^ = 3 we have what can be understood 
as a prototype, (1,1) nanotube model, the one of smallest length L = 1, with periodic 
boundary conditions. For general A^, one can conceive a gedanken device, like the one 
illustrated pictorially for A^ = 6 in Fig. ^d. The A^ vertical lines represent a realization 
of the A sublattice while the B sublattice is mimicked by the central object. The 
radial tracks in the Figure represent conducting paths linking each A site to each B site 
according to the topology of our model. 



z 

^=6 




Figure 6. Physical representation of the CBG for N — 1 (a) and N — 2 (b) which is 
equivalent to a one dimensional ring of length L = 2, 4 respectively. For N = 3 (c) 
we have the (1,1) nanotube of smallest length and periodic boundary conditions. For 
N — 6 (d) we present the gedanken device described in the text. 

Our solution is an example of antiferromagnetic ground state in a model of itinerant 
electrons; it may provide useful hints about antiferromagnetism outside the strong 
coupling regime (where the Hubbard model can be mapped onto the Heisenberg model). 

The model is invariant under an arbitrary permutation of the ^-sites and/or of the 
fi-sites. In the symmetric case |^| = \B\ there is an additional Z2 symmetry because of 
the A ^ B exchange. In what follows we shall focus on the symmetric case and we call 
Af the number of sites in each sublattice, |A| = 2Af. 

The one-body spectrum has three different eigenvalues Eg = —t, £0 = and £g = t 
with degeneracy dg = 1, do = 2J\f — 2 and dg = 1 respectively. We use the convention 
if: > so that Eg is the lowest level and we shall call Shf the set of zero-energy one-body 
eigenstates, \Shf \ = 2Af — 2. The zero-energy one-body orbitals can be visualized by 
a simple argument. Consider any pair i, j, with i 7^ j, of sites belonging to the same 
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sublattice, say A, and a wavefunction ifji^j taking the values 1 and -1 on the pair and 
elsewhere in A and in B. It is evident that ipij belongs to Shf- Operating on ipi^j by the 
permutations of SV we can generate a (non-orthogonal) basis of A/" — 1 eigenfunctions 
vanishing in B\ further, by means of the Z2 symmetry, we obtain the remaining orbitals 
of Shfi which vanish on A. This exercise shows that the group considered above justifies 
the {2N — 2)-fold degeneracy of the one-body spectrum. 

We denote by g {g) the operator which annihilates a particle in the eigenstate with 
energy Eg (eg). Then, the kinetic term K can be written as 

K = -tY^ialoa - gloa) ■ 

a 

Next, we introduce the one-body eigenstate a||0) of Shf with vanishing amplitudes on 
the B sublattice. Similarly, 6||0) has vanishing amplitude on the A sublattice and 
therefore the {2M — 2)-body state 

1$?^) = al . . . aV_i.6L • • • &ki.|0), ^ = (16) 

is an eigenstate of K and of W with vanishing eigenvalue. In the following we shall 
use the wording W = state to denote any eigenstate of H in the kernel of W. It is 
worth noticing that by mapping the ^-sites onto the ;B-sites and viceversa, \^af) retains 
its form except for a spin-flip; we call this property the antiferromagnetic property for 
obvious reasons. 

The state plays a crucial role in building the exact ground state at half 

filling. We first observe that |$af) has non- vanishing projection onto all the total spin 
subspaces S = 0, ...,A/' — 1. Let us denote with |$^^) the singlet component and 
with 1$^^), m = 0, ±1, the triplet component. We further introduce the short-hand 
notation \g'^) = 5'j5'||0) and \g^) = (^jf^jlO) for the two-body singlets that one obtains 
from the lowest and the highest energy orbitals g and g, and m = 0, ±1, for 

the corresponding triplet. Then, one can prove [SHj that the interacting ground state 
|$[/(|A|)) can be written as 

i^uim = [7.1^7°) +7,1/)] ® 1$^^) +70 E i-riigtn ® i^Yf"") , 



b 



where (7^, 7^,70) is the lowest energy eigenvector of a 3 x 3 Hermitean matrix 

We observe that l^^/dAI)) is an eigenstate of the total number operator n°' + n 
of particles in the shell Shf, despite the fact that -|- n*,if] 7^ 0. Such a 
remarkable property {shell population rigidity) relies on the cancellation of those 
scattering amplitudes that do not preserve the number of particles in Shf- We have 
shown [SH] that this cancellation takes place provided in Shf the {2Af — 2)-body state 
is a = state. We have also found that the ground state energy -E(|A|) is negative 
for any value of the repulsion U ; qualitatively, we may say that the particles manage 
to avoid the double occupation very effectively. Furthermore, ii^(|A|) is a monotonically 
increasing function of U and Af due to the existence of non-trivial correlations even 
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for large M . The nature of these correlations has been investigated by computing the 
expectation value of the repulsion. We have shown that for any finite M there is a 
critical value of U yielding maximum repulsion. Even more interesting is the magnetic 
order of the ground state. Due to the SV ® Sj^ ® symmetry, the spin-spin correlation 
function GspmiiJ) = ($t/(|A|)|S'fS'||$[/(|A|)) can be written as 

(Go i= j 

Gspin{i,j) = < Con ieA(ieB) and j eA (j G B) 
[ Gos teAiieB) and j eB (j e A) 
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Figure 7. a) Gon versus U in the range < U < 20 for three different values of the 
number of sites Af — 2, 3, 10. 6) GoS versus U in the range < U < 20 for three 
different values of the number of sites Af — 2, 3, 10. The hopping parameter has been 
chosen to be t = 1 in both cases. 

In Fig. [7|we report the trend of Gon and Gofr versus U for three different values of 
A/" = 2, 3, 10. According to the Shen-Qiu-Tian theorem |59j. Gon is always larger than 
zero while Goaf is always negative. Next, we consider the ground state average of the 
square of the staggered magnetization operator 

mlp ^ -L(^^^l[Y^eit)Sn'\^u), e{{) = 1, -1 hr t e A, B . (17) 
l^^l ieA 

The Shen-Qiu-Tian theorem implies that each term in the expansion of Eq. ()17|1 is 
non-negative. We emphasize, however, that for |^| = \B\ this does not imply that m\p 
is an extensive quantity! Remarkably, in the CBG m\p = Go + (A/" — l)Gon — ■A/'Gofr 
is extensive for any value of the on-site repulsion U and provide the first example of 
antiferromagnetic ground state in a model of itinerant electrons. 

4-2. Half-Filled Hubbard Model on the Square Lattice for f/ = O"*" 
We consider the Hubbard Hamiltonian on the square lattice 

H = K + W=^- J2 (4crV + h.C.) + f/E^rT^ri (18) 
(r,r')<T r 

with r = {ix,iy), ix,iy = 1,...,N\ and the sum on (r, r') is over the pairs of 
nearest neighbors sites. The point symmetry is G4t,; besides, H is invariant under 
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the commutative Group of translations T and hence the space Group G = T ® C^^] 
® means the semidirect product. In terms of the Fourier expanded fermion operators 
(periodic boundary conditions) Ck = ]^ Z^r ""Cr, we have K = Z^k ^kCko-Cko- with 
^k = 2t[coskx + cos ky]. Then, the one-body plane wave state c|^o-|0) = 1^^) 
eigenstate of K. 

In this Section we build the exact ground state of the Hubbard Hamiltonian ()18p 
at half filling and weak coupling for a general even N^- Once a unique non- interacting 
ground state is determined, one can use the non-perturbative canonical transformation 
to test the instability of the system towards pairing; this will be done in the next Section 
for A^A = 4. 

The starting point is the following property of the number operator = cjcr (for 
the moment we omit the spin index). 

Theorem : Let S be an arbitrary set of plane- wave eigenstates {|kj)} of K and 
(nr)jj = (kj|nr|kj) = -i^-e*^''*"'^^)''" the matrix of rij. in S. This matrix has eigenvalues 

Ai = ]^ and A2 = • • • = \\s\ = 0. 

Note that \S\ < N^; if \S\ = the set is complete, like the set of all orbitals, 
and the theorem is trivial (a particle on site r is the rir eigenstate with eigenvalue 1). 
Otherwise, if l^l < N^, the theorem is an immediate consequence of the fact [HI] that 

det|(nr)., - XS.,\ = {-Xf\~\\^ - A), Vr. (19) 

Let Shf denote the set (or shell) of the k wave vectors such that ^k = 0. At half 
filling (iV^ particles) for U = the Shf shell is half occupied, while all k orbitals such 
that ek < are filled. The k vectors of Shf lie on the square having vertices (ivr, 0) and 
(0, ±7r); one readily realizes that the dimension of the set Shf is \Shf \ = 2N\ — 2. Since 
A^A is even and H commutes with the total spin operators, at half filling every ground 

(2 2 \ ^ 
^ _^ j degenerate 

unperturbed ground state configurations with = 0. Most of the degeneracy is 
removed in first-order by W. We shall be able to single out the unique ground state of 
H by exploiting the Lieb theorem. 

The first-order splitting of the degeneracy is obtained by diagonalizing the W matrix 
over the unperturbed basis; like in elementary atomic physics, the filled shells just 
produce a constant shift of all the eigenvalues and for the rest may be ignored in first- 
order calculations. In other terms we consider the truncated Hilbert space Ti spanned 
by the states of Nj^ — 1 holes of each spin in Shf, and we want the exact ground state(s) 
of W in H] by construction TC is in the kernel of K, so the ground state of W is the 
ground state of H as well. Since the lowest eigenvalue of W is zero, it is evident that 
any W = state in 7-^ is a ground state of H. 

To diagonalize the local operator W in closed form we need to set up a local basis 
set of one-body states. If Shf were the complete set of plane- wave states k, the new basis 
would be trivially obtained by a Fourier transformation, but this is not the case. We 
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introduce a set {Iv?^'"-')} of orbitals such that is diagonal in this basis. The eigenvectors 
Iv?^"-*) of nr=o and those Iv^^""-*) of other sites r are connected by translation and also by 
a unitary transformation, or change of basis set. Picking r = e/, / = x means e/ = (1, 0) 
or transfer by one step towards the right and / = y means e/ = (0, 1) or transfer by one 
step upwards. The unitary transformation reads: 



The transfer matrix T/ knows all the translational and point symmetry of the system, 
and will turn out to be very special. 

For large A'a, to find {I'/'a^)} it is convenient to separate the k's of Shf in irreducible 
representations of the space Group G= C^^ ® T. Choosing an arbitrary k G iS/i/ with 
kx > ky > 0, the set of vectors -Rjk, where Ri G C^v, is a (translationally invariant) basis 
for an irrep of G. The high symmetry vectors (0, vr) and (tt, 0) always trasform among 
themselves and are the basis of the only two-dimensional irrep of G, which exists for any 
A^A- If is even, one also finds the high symmetry wavevectors k = (±7r/2, ±7r/2) 

which mix among themselves and yield a four-dimensional irrep. In general, the vectors 
Rih are all different, so all the other irreps of G have dimension 8, the number of 
operations of the point Group C^^. 

Next, we show how to build our local basis set and derive W = states for each 
kind of irreps of G. For illustration, we shall first consider the case A^a = 4; then 
Shf contains the bases of two irreps of G, of dimensions 2 and 4. The one with basis 
k^ = (vr, 0), k^ = (0, tt) breaks into Ai © Bi in 

The eigenstates of (?T,r=o)ij = (kj|rar=o|kj), with i,j = A,B, are = + 

Ik^)) with Ai = 1/8 and \iP'b^) = -^H^a) — \^b)) with A2 = 0. Since under translation 
by a lattice step Ti along the I = x,y direction |k) — > e^'^'lk), using Eq. (pUj) one finds 
that {i/Ja-^) ^ {—ly'i \ip'^^) , with 9'^ = 1, 9y = 0; so \iI>a-^) has vanishing amplitude on 
a sublattice and IV'^^) on the other. The two-body state |V'Ai)o-|V'bi)-o- has occupation 
for spin a but not for spin —a on the site r = 0; under a lattice step translation it flips 
the spin and picks up a (-1) phase factor: |V'Ai)o-I^Si)-cr ^ Wbi) ^Wai) -o-'^ therefore it 
has double occupation nowhere and is a 14^ = state (W = pairj^^.j^). 

The 4-dimensional irrep with basis ki = (7r/2,7r/2), k2 = (— 7r/2, 7r/2), ks = 
(7r/2, -7r/2) k4 = {-it/2, -n/2) breaks into Ai ® B2 ® E in C^„\ letting / = 1,2,3,4 
for the irreps Ai, i?2, E^, Ey respectively, we can write down all the eigenvectors of 
{nr=o)ij = {ki\nr=o\kj) , with i,j = 1, ... ,4, as {ipj) = Yli=i O/Jki), where O' is a 4x4 
orthogonal matrix. The state with non-vanishing eigenvalue is again the Ai eigenstate. 
After a little bit of algebra we have shown jBlj that under Ti the subspace of Ai and 
B2 symmetry is exchanged with the one of E^ and Ey symmetry. Thus we can build 
a 4-body eigenstate of W with vanishing eigenvalue: \4''Ai4''B2)<^\'^'Ej:i^'Ey)-(T- before 
under a lattice step translation this state does not change its spatial distribution but 



a ^ -a without any phase factor: WA,fB2)'^\^k^'Ey)-^ ^ lV^k^E^)<^l^^iV'B2)-a- 




(20) 
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Now we use these results to diagonahze nr=o on the whole set Shf (we could have 
done that directly by diagonalizing 6x6 matrices but we wanted to show the general 
method). The eigenstate of nr=o with nonvanishing eigenvalue always belongs to Ai. 
The matrix rip has eigenvalues 3/8 and (5 times) 0, as predicted by Eq. ()19|). For 
r = the eigenvector of occupation 3/8 is Iv^i^'*) = -^Wa^) + \I\Wa^)- The other 
Ai eigenstate of nr=o has eigenvalue and reads: |v52°^) = \[\Wa^ ~ '^Wa^)- The 
other eigenvectors, whose symmetry differs from Ai, are \^^^^) = iV'sa)' ~ Wbi)^ 

lv'5°'') = WeJ) Iv^e^"*) = Wsy)- finds |UT| that the transfer matrices Ti of Eq. 
such that = J2j are antiblock diagonal. Thus, the local basis at any 

site r splits into the subsets Sa = U'^i'^), l^^f^), ly^f^)}, and Sb = {Iv'i^''), I'^f^), Iv^e^"*)}; 
a shift by a lattice step sends members of Sa into linear combinations of the members 
of Sb, and conversely. 

Consider the 6-body eigenstate of K 

In this state there is partial occupation of site r = with spin a, but no double 
occupation. It turns out that a shift by a lattice step produces the transformation 

\^AF)a < > -\^AF)-a 

that is, a lattice step is equivalent to a spin flip, a feature that we have already met 
in Section ITTl ( antiferromaanetic property). Since the spin-flipped state is also free of 
double occupation, |$af)o- is a = eigenstate. A ground state which is a single 
determinant is a quite unusual property for an interacting model like this. 

Note that |v5i°V2°^) is equivalent to I^AiV'Ai)' because this is just a unitary 
transformation of the Ai wave functions; so |$AF)cr can also be written in terms of 
the old local orbitals (without any mix of the local states of different irreps of G): 

I^Af). = \^AA^'B,).\^'k^'Ej'Ej--- (21) 

This form of the ground state lends itself to be generalised to arbitrary even A^a, see 

Refs.[nn],inii. 

A few further remarks about I^af)^ are in order. (1) Introducing the projection 
operator Ps on the spin S subspace, one finds that Ps\^AF)a = \^af)(^ = 
0, . . . , A^A - 1. Then, ^{^AF\W\^AF)a = Es=I^ A^af\W\^^af)^ = 0, and this implies 
that there is at least one ground state of in ?i for each S. The actual ground state 
of H at weak coupling is the singlet I^'af)'^- (2) The existence of this singlet W = 
ground state is a direct consequence of the Lieb theorem |25j. Indeed the maximum spin 
state \^AF~^)cr is trivially in the kernel of W; since the ground state must be a singlet 
it should be an eigenvector of W with vanishing eigenvalue. (3) The above results and 
Lieb's theorem imply that higher order effects split the ground state multiplet of H and 
the singlet is lowest. (4) The Lieb theorem makes no assumptions concerning the lattice 
structure; adding the ingredient of the G symmetry we are able to explicitly display the 
wave function at weak coupling. 
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Energy 


Irrep of Q 


Degeneracy 


4 




1 


2 


A4 


4 





Q4 


6 


-2 


Ai 


4 


-4 




1 



Table 4. One-body spectrum of the 4x4 Hubbard model for t — —1. We have used 
the notation introduced in Ref. j(i5j in labelhng the irreducible representations. 

Using the explicit form of ^5=0 one finds that Ps=o\^AF)a = —Ps=o\^AF)-a- This 
identity allows us to study how the singlet component transforms under translations, 
reflections and rotations. In particular the antiferromagnetic property tells us that the 
total momentum is Ktot = (0, 0). To make contact with Ref. [52] we have also determined 
how l^^i?) transforms under the C^^ operations with respect to the center of an arbitrary 
plaquette. It turns ouffH!] that it is even under reflections and transforms as an s wave 
if N\/2 is even and as a d wave if N\/2 is odd. 

In the next Section we use these results, togheter with the non-perturbative 
canonical transformation, to study the doped 4x4 lattice at half filling. Since the non- 
interacting ground state at half filling is now well known and unambiguously defined, 
the expansion ()12|) can be performed in a unique way. 

4-3. Pairing in the Doped Hubbard Antiferromagnet 

The one-body spectrum of the 4x4 Hubbard model has 5 equally spaced levels, see 
Table |3 The space Group G [containing the translations and the 8 (74^, operations] can 
not explain the degeneracy 6, since in the 4x4 lattice the largest dimension of the irreps 
is 4. As observed by previous authors [HSljinil, the 4x4 lattice can be mapped into the 
2x2x2x2 hypercube since each pair of next to nearest neighbor sites has two nearest 
neighbor sites in common. This implies that H is invariant under a new and largest 
symmetry Group; let us call it Q. 

Due to the importance of the symmetry in our configuration interaction mechanism, 
we have explicitly calculated the character table of Q taking into account an extra non- 
isometric symmetry operation [HHl- Q has 384 elements, 20 classes and hence 20 irreps 
whose dimensionality fully justifies the degeneracies of Table IH 

As observed in Section 14.11 the canonical transformation applies when two holes are 
added to a non-degenerate vacuum. To study the system at and close to half filling, we 
have to use the results of Section 14.21 In the following we will solve the problem of two 
electrons added to the half filled system. 

4-3.1. W = Pairs In order to study the W = pairing in the doped 4x4 
antiferromagnet, we consider W = pairs in the six-fold degenerate one-body level 
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[belonging to Vt^. Exploiting the W = Q theorem, we have found[nH] W = Q pairs with 
symmetry Fi, S2, f^i and triplet pairs with symmetry S3, 0.2, VL^. Here, we are 
using the notation of Ref. |65j: the irreps VL have dimension 6, the S's have dimension 3 
and r's have dimension 2. 

Exact diagonalization results[32., show that for U/t < 3 and 16 — 2 = 14 holes 
the ground state is sixfold degenerate. Below, we use the canonical transformation and 
prove that the ground state corresponds to an fii = electron pair over the half-filled 
system. For U/t > 3 and the same number of holes a level crossing takes place: the 
ground state is threefold degenerate and contains a state with momentum (0, 0) and a 
doublet with momentum (vr, 0) and (0,7r). Again, the computed ground state can be 
assigned to a S2 electron pair over the half filled svstemj65j. 

4-3.2. Pairing mechanism We consider the ground state of the 4x4 Hubbard model 
with 14 holes; aside from the 10 holes in the inner Ai and Ai shells (see Table Ej), the 
outer 0,4^ shell contains 4 holes. We intend to show how pairing between two electrons 
added to the antiferromagnetic 16-holes ground state (half filling) comes out. We use 
the antiferromagnetic ground state |$^^°) as the non- interacting ground state of the 
configuration interaction expansion (fT^ . With respect to this electron vacuum, the m 
states are now W = pairs of electrons added to |$;^^°). In the 4x4 lattice the one-body 
energy levels are widely separated and the dominant interaction is between electrons in 
the same shell. Therefore, we consider as m states only the W = electron pairs in 
the shell Shf, i.e., those belonging to the irreps fli and S2, and we neglect the high- 
lying unoccupied orbitals. Explicit calculations [021 show that the effective interaction is 
attractive for both Qi and S2 = electron pairs and that the corresponding binding 
energy is A^^ = —61.9 meV and A^j = —60.7 meV for f/ = — t = 1 eV. Therefore, the 
weak coupling ground state can be interpreted as an Qi W = electron pair over the 
antiferromagnetic ground state. This result agrees with exact diagonalization data|32j. 

5. Carbon Nanotubes and Triangular Cobalt Oxides 
5.1. Nanotubes 

There is experimental evidence that the critical temperature Tc in alkali-graphite 
intercalation compounds C^M (where M is a given alkali metal) grows as x decreases j66j. 
Under high-pressure, high metal concentration samples such as CgK, C3K, C4Na, CsNa, 
C2Na, C2Li have been synthesized; for C2Na the value of is 5 K while for C2Li, 
Tc=1.9 K. Recently, Potassium|67j and Lithium [H^ have been intercalated also in 
ropes of single- and multi-wall carbon nanotubes (the highest metal concentration was 
obtained with Lithium in C2Li), and a net charge transfer between the alkali- metals 
and the carbon atoms has been predicted [Snj- The alkali-metals cause little structural 
deformation, but increase the filling of the original bands. Nanotubes close to half 
filling are deemed to be Luttinger liquids down to milli-Kelvin temperatures |7()j .171 j. 
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Here|72j, we use the Hubbard Hamiltonian H on the wrapped honeycomb lattice to 
represent the valence bands of single-wall carbon nanotubes (SWNT) and we apply our 
symmetry- driven configuration interaction pairing mechanism based on the existence 
of = pairs. We shall focus on armchair (A^, N) SWNT and show that the pair 
binding energy grows as the number of electrons per C atom increases. Furthermore, 
a stronger binding in nanotubes than in graphite sheets is predicted which suggests a 
higher critical temperature in the former. This is also supported by the measurements 
of a Tc ~ 15 K in the 4 Angstrom SWNT by Tang et al.^. 

Any armchair SWNT has a bonding (-) and an antibonding (+) band, and the 
Fermi line has C2v symmetry (Ce?, is the symmetry group of the graphite sheet). Since 
these systems are usually doped with electrons, we take the Fermi level ep io lie 
in the antibonding band. We use the W = Q theorem for electrons with opposite 
quasimomentum and we find W = Q pairs belonging to the pseudoscalar irrep of 
Let (a, 6) denote the basis of the Bravais lattice and m (k, ^) the periodic part of 
the Block function of quasi-momentum k, with ( = a,b. The singlet pair wavefunction 
reads [72] 

(k, Ri, R2) = sin {K{X, - X2)) X 

X u* (k, Ci) u* (-k, C2) e^^''^^-^^) - u* (k, C2) u* (-k, Ci) e-'^^^^-^^); , 

with Rj = (Xj, Yj) the origin of the cell where the particle i lies. We can verify by direct 
inspection that the W = pair wavefunction ip^^^^ vanishes for Xi = X2, that is if the 
particles lie on the same annulus of the armchair tube. 

The effective interaction Wes between the particles of aW = pair can be obtained 
analytically by the canonical transformation approach described in Section We let 
e(k) be the one-body energy excitation with momentum k in the antibonding band and 
we call P/4 a quarter of the empty part of the FBZ. The effective Schrodinger equation 
for the pair reads 

[2e{k) + WF + F{k,E)]a^+ J2 W^s {Kk' , E) = Ea^ , (22) 

k'eI'/4 

where Wp is the first-order self-energy shift (which we found to be independent of k) and 
F(k, E) is the forward scattering term (which does not contains any direct interaction 
between the particles of the pair). Eq. requires a self-consistent calculation of E 
(since W^s and F are i?-dependent). We show below that E = 2e p + Wf + F^i^(kp) + A, 
with a positive binding energy —A of the W = pair; here Finin(kir) is the minimum 
value of -F(k, E) among the k^- wave vectors on the Fermi line. 

First, we got a direct verification that pairing actually occurs in the (1,1) nanotube 
of length L = 2 (in units of the lattice spacing) and periodic boundary conditions. As 
for the CUO4, we compute the quantity A(4) = -E'(4) + E{2) — 2E{3) and we compare 
it with A (obtained from the canonical transformation). A (4) has been computed in a 
large range of U/t values, and its trend is shown in Fig. |H^a). In Fig. IH^b) it is reported 
the comparison between A (4) and A. We observe that the analytical value |A| is ~ 2 
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Figure 8. (a) Trend of A(4) in units of t versus Log U/t in the range -1-^3 for the 
(1,1) nanotube. (b) Comparison between A and A(4). 

times greater than |A(4)| for f//t ^ 1, as for the CUO4 cluster. However, the analytical 
approach predicts the right trend of the binding energy and it is very reliable in the 
weak coupling regime. 

Next, we consider supercells of 2 N x L = Nc cells, where L is the length of the 
(A^, A^) nanotube in units of the lattice spacing. We solve the Cooper-like equation 

in a virtually exact way for up to 6 and L up to 25, using U/t = 2.5 (which is 
of the correct order of magnitude for graphite [Z3],[Z1])- The canonical transformation 
overestimates A in this range of U/t, but remains qualitatively correct. The calculations 
are performed with the Fermi energy ep varying between 0.8 t and 1.1 t (half filling 
corresponds to ep = 0). As in the (1, 1) nanotube, the W = pairs get bound once 
dressed by the electron-hole excitations. However, the binding energy —A decreases 
monotonically both with the radius and the length of the tube. 

With supercell sizes Nc > 300 numerical calculations become hard and the AEI 
scheme is used in order to get reliable extrapolations. The AEI Ves remains fairly stable 
around ^ 1.5-7-2 t for N > 2 with increasing L. Furthermore, V^g is largely independent 
of the Fermi energy. The weak dependence of Ves on the lenght L allows for extrapolating 
the asymptotic value of the binding energy Aasympt(A^) = liniL^oo A(A^, L), see also 
Section ing The results are shown in Fig. Efa) with V^fj = 1.5 t. We found that Aasympt 
is strongly dependent on the filling at fixed A^; the sharp maximum at the optimal 
doping Sp ^ t (which corresponds to a number of electrons per graphite atom of 1.25) 
can be understood in terms of a corresponding peak in the density of states. In the 
optimally doped case — Aasympt (^) decreases monotonically as the radius of the tube 
increases, see Fig. El^b). The decreasing of the binding energy with A^ is suggested by 
recent measurements on nanotubes with diameter of few AngstromjS]. However, in the 
limit of large A^, Aasympt (^) remains stable around 0.0028 t and may be interpreted as 
the binding energy of the = pair in an optimally doped graphite sheet. 

We emphasize that our pairing mechanism uses degenerate electronic states that 
exist in 2d (or quasi 2d) and works away from half filling. A proper account for the 
transverse direction is crucial in order to have a non-Abelian symmetry group and hence 
W = Q pairs. The ^/'^'^^l p^ir wavefunction vanishes when the transverse component 
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Figure 9. (a) Results of the canonical transformation approach with U/t ~ 2.5. 
— Aasympt as a function of the Fermi energy ep for iV = 4 (black boxes), = 6 (empty 
triangles) and iV = 10 (grey diamonds). The Fermi energy varies in the range 0.8-^ 1.1 
t. (b) —Aasympt as a function of N for N in the range 6-=-36 with ep = t and Average 
Effective Interaction V = 1.5 t. In both figures — Aasympt is in units of t. 



ky = 0. This opens up the interesting possibihty that in nanotubes two distinct 
superconducting order parameters appear in the phase diagram, if it turns out that 
close to half-filhng there is another one which breaks down the Luttinger hquid. 

5.2. Triangular Lattice in W = Theory 

We have also considered a symmetric triangular Hubbard lattice, which may be 
relevant to the newly discoveredjS] Nai.Co02 superconductors which are now exciting 
considerable interest [ZHl , [ZH! • The 7 atom centered hexagonal cluster with open 
boundary conditions yields no pairing for any filling. When opposite sites of this 
cluster are identified, one obtains a 4-site cluster which is the smallest one with periodic 
boundary conditions. With a hopping integral t = — 1 we find that A (4) is negative and 
shows a similar trend versus U as in the CUO4 case; the pairing energy exceedes O.lt for 
U ~ 4|t|. 

6. Superconducting Flux Quantization 

Bulk superconductors quantize the flux through a hole in half-integer multiples of the 
fiuxon 00, because the quasiparticles that screen the vector potential carry charge 2e. 
In finite systems the signature of superconductivity is a ground state energy minimum 
at = that is separated by a barrier from a second minimum at = 0o/2- With 
increasing the size of the system, the energy (or free energy, at finite temperature) barrier 
separating the two minima becomes macroscopically large, and bulk superconductors 
can swallow up only half integer numbers of flux quanta. As emphasized by Canright and 
Girvin^j, the flux dependence of the ground state energy is definitely one of the most 
compelling ways of testing for superconductivity, and the existence of the two minima 
separated by a barrier is a strong indication of superconducting flux quantization. 
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6.1. General Group Theory Aspects of Gu-0 Systems 

In the present problem, with a repulsive Hubbard model, the mechanism of attraction 
is driven by the C^^ symmetry, and cannot operate in an unsymmetric geometry. The 
flux must be inserted in such a way that the system is not distorted. In the following we 
consider Cu-0 systems with C^^ symmetry with respect to a central Copper ion, and 
insert the magnetic flux in such a way that only 4 central triangular plaquettes feel a 
magnetic field (see Figs. El and ITT|l and the rest of the plane only experiences a vector 
potential. 

Consider the pattern of Fig. [TUl Here, the Cu sites are marked by black dots and the 
Oxygen sites by empty dots; the X's stand for tubes carrying flux each, symmetrically 
disposed around the central Cu. According with the Peierls prescription we include the 

o o o o o 
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Figure 10. Pattern of the veetor potential A due to 4 flux tubes (crosses) carrying 
flux (f). black dots stand for Cu atoms and empty dots stand for O sites. The line 
integral of A along each bond parallel to the arrow is ^ . 



magnetic effects in the Hubbard model by setting 

where 0o = 27r/e is the flux quantum and j and f are the position of the two lattice 
sites connected by tjji. Varying (j) by an integer multiple of 0o corresponds to a gauge 
transformation leaving all the physical properties unchanged. The arrows help to 
visualise a convenient choice of the gauge at general (f). Namely, running along an 
oriented bond in the sense of the arrow, A ■ dr = |; along the other Cu-0 bonds, 
not marked in the figure, / A ■ = 0. One sees that in this way the flux through any 
closed path corresponds to the number of tubes surrounded by the path. The reflection 
operations of C^^ are equivalent to ^ —0, reverse the directions of the arrows and 
for a generic the symmetry Group reduces to ^4. However, at = ^ the reversal of 
the magnetic field in the tubes corresponds to a jump by 0o, and this is equivalent 
to a gauge transformation: this implies that the symmetry Group gets larger, the 
new symmetry operations being reflections supplemented by a gauge transformation. 
Indeed, the hopping parameter tjji becomes itjji along the arrows, while it remains 
equal to tjji along the unmarked bonds of Fig. (TnT a). Any reflection operation simply 
changes the signs of all the hoppings along the marked bonds. Now consider the unitary 
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transformation S which changes the signs of all the Cu orbitals along both diagonal, 
except the central Cu. Since S also has the effect of reversing all the arrows, a x is 
a symmetry, for all reflections a in (74^,. Moreover, since the product of two reflections 
is a rotation, the Group including the rotations and the reflections multiplied by 
S is isomorphic to C4„. The W = Q pairs appropriate for half a flux quantum must 
involve two holes belonging to the degenerate irrep of C*4„. In this way, at = ^ the 
full symmetry is restored, allowing again for pairing and negative A. If pairing also 
occurs at = the superconducting flux quantization arises from a level crossing 
between the ground state associated with the paired state at = and the ground 
state associated with the paired state at = 

6.2. Application to the CuO^ Case and Numerical Results 

As an illustrative application of the previous symmetry argument, let us investigate one 
more time the CUO4 cluster. We expect that this system is a very good candidate to 
exhibit superconducting flux quantization, since it hosts two W = Q pairs of different 
symmetries. As discussed above, this condition is necessary for the development of a 
level crossing. 
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Figure 11. (a) Ground state energy of the CUO4 with four holes as a function of 0/0o- 
For (t)/(t>o between and 1/4 and between 3/4 and 1 the ground state has B2 symmetry; 
for (l)/<f)o between 1/4 and 3/4 it belong to Ax. Here t — leV, tpp — — O.OleV, U — 5eV; 
i?(4, (/)) is in eV. (b) Topology of the CUO4 cluster in presence of 0; the crosses stand 
for flux tubes. 



In agreement with the previous prescription, we have to insert 4 flux tubes around 
the central Copper [see Fig. IHfb)] in order to introduce a closed path around the 
center, where screening currents can respond. Every 0-0 bond collects the Peierls phase 
^ / A ■ (ir = 2iTi-^; by symmetry, t is unaffected by the flux. Thus, the Hamiltonian 
reads as in Eq. ^ with 

We observe that a flux of the order of a fluxon in a macroscopic system would be a 
small perturbation; in the small cluster, however, the perturbation is small only if the 
hopping integral \tpp\ is taken small compared to t. Numerically, the computations were 
performed with tpp = —0.01 eV. By exact diagonalization, we have found that the ground 
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state energy E{(j)) of the CUO4 with 4 holes, as a function of 0, has clearly separated 
minima at zero and half a flux quantum [see Fig lTlT a)]. Furthermore, A (4) is negative 
both at = and = 0o/2: A(0 = 0) = -43meV and A(0 = 0o/2) = -32meV. 

The physical interpretation of Fig. ITlT a) is the following. When the magnetic 
flux is inserted into the system, the W = Q pair of B2 symmetry creates a diamagnetic 
supercurrent that screens the external field. Such a current flows through the 0-0 bonds 
and form closed loops. As grows the energy of the system also increases, signaling that 
the W = Q pair is spending its binding energy to screen the field. At a quarter fiuxon a 
level crossing occurs, producing a second minimum at = 0o/2. Here the Hamiltonian 
is real again and the Ai pair is energetically favoured. As the flux increases further, the 
Ai pair produce a new diamagnetic supercurrent until the initial situation is restored at 
= 00. The pairing state at zero flux and half fiuxon are orthogonal. There is a clear 
analogy with the BCS theory; in that case, the Cooper wavefunction has s symmetry 
and the total magnetic quantum number of the pair vanishes in the absence of flux, but 
not at half a flux quantum. 

It is worth to note that if the Hubbard repulsion U is set to zero, the second 
minimum at half a fiuxon disappears and a trivial paramagnetic behaviour is observed. 
This is a further evidence that the superconducting behaviour of the system is induced 
by repulsion. 

6.3. Rings of Symmetric Clusters 

We have also studiedjTS], [7^1,1101 bound pair hopping and Superconducting Flux 
Quantization (SFQ) in systems with CUO4 units as nodes of a graph A. For such 
systems the Hamiltonian is i^tot = Hq + Hr with 

while Ht- is an intercell hopping Hamiltonian. Here, is the creation operator for a 
hole onto the Oxygen i = 1, .., 4 of the a-th cell and so on. The point symmetry group of 
Hq includes S]f^\ with |A| the number of nodes. We shall consider an intercell hopping 
which preserves the S4 subgroup of in order to have W = pair solutions. 

6.3.1. 0-0 Hopping and SFQ Consider a hopping term that allows a hole in the i-th 
Oxygen site of the a-th unit to move towards the i-th Oxygen site of the /3-th unit with 
hopping integral r^^ = Itq^Ic*^"'': 

= E E Pi,iaPl3,ia ■ 

a,/3eA icr 

For = 2|A| and = 0, the unique ground state consists of 2 holes in each CUO4 
unit. The intercell hopping produced by small Itq/jI <^ |A(4)| allows for studying the 
propagation of p pairs added to a background 2|A| holes. When U/t is such that 
A (4) < 0, each pair prefers to lie on a single CUO4 and for A^ = 2|A|-f2p the unperturbed 
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ground state is 2^x(^^^ times degenerate (since each CUO4 can host two degenerate 
W = Q pairs). In the low-energy singlet sector, the problem is solved analytically to 
second order in and mapped into an effective Hamiltonian for p hard-core bosons 
with a complex effective hopping integral J that we have calculated analitically and 
studied as a function of the ratio U/t. For ring-shaped systems, the effective model is 
equivalent to the Heisenberg-Ising spin chain governed by the Hamiltonian 

^Hi = E J[2w>a+i + + e-wia>-+J (23) 

a=l 

where the cr's are Pauli matrices, spin up represents an empty site and spin down 
represents a pair, r] is the so called anisotropy parameter and in our case rj = —1. For 
= 1, we have the isotropic Heisenberg interaction. By performing a Jordan- Wigner 
transformation, the Hamiltonian in Eq. (j23|) can also be mapped into a model of spinless 
fermions on the ring. In the absence of a threading magnetic field {(f) = 0) the problem 
was originally studied by Bloch|^3 and then exactly solved by Hulthen|H2] [in the case 
T] = —1] and Orbach|HHl [in the case r] < —1] using the Bethe's hypothesis |H1]. A 
systematic analysis in the whole range of parameters was given by Yang and Yang in 
a self-contained series of papers|SS]- Here we just recall that the model has a gapless 
phase if \r]\ < 1, corresponding to the conducting state, while an insulating phase sets 
in for ?7 < — 1. As in the Id Hubbard model, the "magnetic perturbation" (0 7^ 0) 
does not spoil the integrability and the Heisenberg-Ising Hamiltonian remains exactly 
solvable by the Bethe-ansatz method. Let us write an eigenfunction of Hm as 



p 



where P is a permutation of the integers 1, . . . ,p and Ap are p\ coefficients. Shastry 
and Sutherland [HS] have shown that the variables kj are given by 

\A\kj = 2nlj + 47r^ - ^ 0{kj, h) (24) 
with a phase shift 



d{k,q) = 2tan"^ 



rism[{k — q)/2] 



cos[{k + q) /2] - r] cos[{k - q) /2] 

From Eqs. (j2Sl)-(|211) we readly see that the ground state energy of the low-energy 
effective Hamiltonian ifni is periodic with period 0o/2, independent of the number of 
added pairs. Thus we conclude that the purely repulsive Cu04-Hubbard ring threaded 
by a magnetic field quantizes the flux in a superconducting fashion if the number of 
particles is 2|A| + 2p with < p < |A|. 

For rings of 2 and 3 Cu04-sites, we also confirmed the analytic results by exact 
diagonalization; for the three-Cu04 ring we observed the SFQ with total number of 
holes 2|A| -|- 2 = 8. The switching on of the hopping r between the Oxygen sites breaks 
the symmetry group Cs^ ® Sf into C^^ ® 'S'4 for real r; in a magnetic field (complex 
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r), this further breaks into C3 (g> 5*4. Real r hfts the degeneracy between the = 
subspace and the subspaces k = 1 and = 2 of C3 (as usual k is related to the crystal 
momentum p = 2TThk/3 in this case), but cannot split k = 1 and 2 because they belong 
to the degenerate irrep of C^^; complex r resolves this degeneracy. 

The three-Cu04 ring is the smallest ring where we can insert a magnetic flux by 
T = |r|e*^, 9 = ^(0/0o)- The energies of the three ground-state multiplet components 
are reported in Fig. IT^ a) for |r| ^ 1^(4) | and U = 5t. At = the ground state 



-8.9882 



-8.9883 




0.3 



-0.3 

















■ 






'I 
■ 




fk=0 


\ 1 


k=2 




\ 1 
\l 


k=l 




■X 
■ N 






l\ 
1 ^ 






■ 

■ J 






■ . 





























0.5 



(a) (b) 

Figure 12. (a) Numerical results for the low energy states of the three-Cu04 ring, as 
a function of the concatenated magnetic flux. The energy is in units of t. (b) Total 
current for the three-Cu04 ring, as a function of the magnetic flux. The current is in 
units oie\T\/h. In both figures, U ^5t, [A(4) w -0.04258 1] |t| = 0.001 



belongs to the k = subspace, while the first excited levels have k = 1 and 2. Their 
spatial degeneration is fully lifted: the k = 1 level increases while the k = 2 level 
decreases up to = 0o/2. As increases, the ground state energy grows quadratically 
in (diamagnetic behaviour). Near = 0o/4 we find a level crossing between k = 
and k = 2, while at0 = 0o/2,A; = O becomes degenerate with k = 1 and the ground 
state energy is in a new minimum belonging to the k = 2 subspace: a sort of "restoring" 
of the = situation is taking place as in the BCS theory [HZI. Indeed, at = 0o/2 the 
symmetry group is C^y 5*4 where C^y is isomorphous to C^^ (reflections a are replaced 
by ag, where g is a suitable gauge transformation). 

Fulfilling the conditions A (4) < and |r| ^ |A(4)|, we have varied U and |r| 
and found analogous trends for the ground state energy. Increasing |r| with fixed A(4) 
lowers the central minumum and depresses the two maxima. On the other hand, if 
|A(4)| decreases at fixed |r| the central minimum and the side peaks are affected in a 
similar way. This is reasonable since the perturbative parameter is |r|/|A(4)|. 

The current operator / = yieldsjHH] a gauge invariant average I and in Fig. 
IT^ b) we display / as a function of in the three ground state sectors k = 1, 2 and 3. 
The current is proportional to the flux derivative of the ground-state energy [see Fig. 
IT^ a)] according to the Hellmann-Feynman theorem and clearly shows a superconducting 
behaviour. 

The ground state energy in each k sector for the non-interacting {U = 0) CU3O12 
ring shows no pairing in CUO4 and indeed the ground state energy is linear in the field at 
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small fields (normal Zeeman effect). The lowest state is = 2 throughout. The absence 
of SFQ is a further evidence of the repulsion-driven pairing mechanism discussed above. 
Thus, the dressed W = Q pair screens the vector potential as a particle with an effective 
charge e* = 2e does. At both minima of -E'(0) we have computed for the CU3O12 ring 
A(8) ^ -IQ-H. 

6.3.2. Cu-Cu Hopping: infinite effective mass The three-Cu04 ring with 2|A| +2 = 8 
holes contains a,W = Q pair which in the case A (4) < gets bound; an intercell hopping 
Tcu is now assumed between Cu sites only. This is perfectly able to carry a one-body 
current when the full system is threaded by the flux, and it does in the absence of 
interactions. However the bound pairs have an interesting behaviour. |rcu| ^ l^(4)| 
produces very much smaller effects on the energy than r does, but anyhow the system 
behaves as a paramagnet and there is no flux-induced splitting of the three k levels. This 
happens because the W = Q pair is strictly localized by the local symmetry. Indeed the 
S'4 label of each CUO4 unit is a good quantum number. No SFQ is observed because the 
screening of the magnetic fleld by the bound pair is forbidden. The superconducting 
pair behaves as if it had a very large effective mass. 

7. W = Q Pairing and Electron-Phonon Interactions 

The role of electron-phonon (EP) interactions in determining the superconducting 
correlations in the Cu-0 planes of cuprates is a very controversial issue. Possibly, 
the pairing mechanism has a predominantly electronic origin, but many high-Tc 
compounds exhibit a quite noticeable doping-dependent isotope effect, suggesting that 
EP interactions could be important and should be included in the theory. In particular 
there is experimental evidence that the half breathing Cu-0 bond stretching mode at 
k = (vr, 0), (0, tt) is signiflcantly coupled with the doped holes in the superconducting 
regime and its contribution may be relevant for the dx2-y2 pairing [5^. 

Here we further investigate this issue by addressing the question if the W = 
pairing available in the Hubbard model survives when the lattice degrees of freedom 
are switched on. When lattice effects are introduced in this scenario, several questions 
arise. In the conventional mechanism, phonons overscreen the electron repulsion; what 
happens if electronic screening already leads to pairing? It is not obvious that the 
phonons will reinforce the attraction while preserving the symmetry. More generally, 
some vibrations could be pairing and others pair-breaking. To address these problems we 
use an extension of the Hubbard model in which bond stretchings dictate the couplings to 
the normal modes of the C4t,-symmetric conflguration, generating a long-range (Frohlich) 
EP interaction. 

We start from the Hubbard model with on-site interaction U and expand the 
hopping integrals tij{ri,rj) in powers of the displacements pi around a C4^-symmetric 
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equilibrium configuration 

a 

where a = x,y. Below, we write down the in terms of the normal modes g^i^: 
P? = J2r)iy S'^vi}) where rj is the label of an irreducible representation (irrep) of the 
symmetry group of the undistorted system and u is a phonon branch. 

Thus, treating the Cu atoms as fixed, for simplicity, one can justify an electron- 
lattice Hamiltonian: 

Hel-latt = Hq + Vtot ■ 

Here Hq = Hq + Hq accounts for the kinetic part of the electron-phonon system and it 
is given by 

the a;'s are normal mode frequencies. The interacting part Vtot = V + W contains the 
Hubbard repulsion W and the electron-phonon interaction V. The latter can be written 
in terms of the parameters $,ri,u = \y\J 2h^~-' where the A's are numbers of order unity 
that modulate the EP coupling strength, while M is the Oxygen mass. Then, 

and the H^f^u operators are given by 

This is physically more detailed than the Hubbard-Holstein model, where electrons are 
coupled to a local Einstein phonon and the superconducting phase has been investigated 
in detail|nn],|ni]- Indeed, in the present context the Hubbard-Holstein model is not fully 
satisfactory because it restricts to on-site EP couplings, which would be impaired by a 
strong Hubbard repulsion. On the other hand the use of a Frohlich-like phonons was 
suggested |n21,|nS],|ni| for modeling the Cu-0 planes in the strong EP coupling regime: 
a long-range EP coupling removes the problem of polaron self-trapping, otherwise 
present in the case of the Holstein interaction, where an unphysically large polaron 
(and bipolaron) mass occurs. 

By generalizing to H^i-iatt the canonical transformation proposed for the electronic 
part H = Hq + W, one can derive an effective interaction between the particles 
in the pair. We obtained a new Cooper-like equation impair] 92) = E\(p) with an effective 
two-body Hamiltonian, acting upon the dressed W = pair \ip). As in Section 3 the 
pairing criterion involves the properly renormalized Fermi energy Ep . if the lowest 

( R) 

energy eigenvalue E is such that E = Ep + A with negative A, the dressed W = pair 
gets bound in the many-body interacting problem and the system undergoes a Cooper 



dtijirj.Vj) 



dtijirj.Vj) 
drf 



dt 



\ iciaCj,a + h.C. 

n J 
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instability. We observed that this extended approach is very accurate at weak couphng 
and is quahtatively predictive also in the intermediate coupling regime. 

As an illustrative application we focussed again on CUO4. This cluster represents 
a good test of the interplay between electronic and phononic pairing mechanisms since 
we can compare exact diagonalization results with the analytic approximations of the 
canonical transformation. CUO4 allows only the coupling to phonons at the centre or at 
the edge of the Brillouin Zone; however, phonons near the edge are precisely the most 
relevant onesjEH]- We recall (see Section 1221) ^^at the pure Hubbard CUO4 cluster with 
0-0 hopping tpp = yields A (4) < 0, and a degenerate W = bound pair with Ai and 
B2 components. A standard Jahn- Teller calculation in which the degenerate ground 
state wave functions interact with the vibrations predicts distortions that already at 
moderate EP coupling destroy the pairingf96]. We also explored the scenario beyond 
this approximation. The analytical solution of the Cooper-like equation shows that the 
electronic pairing is enhanced by the phonon contribution. In particular the binding 
energy of the pair in the B2 symmetry channel increases more rapidly than in the Ai 
channel. 

In order to go beyond the weak-coupling regime we have exactly diagonalized the 
Hamiltonian in a truncated Hilbert space. We obtained several indications. First, 
the catastrophic distortions predicted by the Jahn- Teller approximation are not borne 
out by more complete approaches involving a broader spectrum of electronic states. 
Pairing prevails also at strong coupling in part of the parameter space, in the symmetry 
channels where W = pairs occur. The correct trend is predicted by the canonical 
transformation approach, which also explains the pairing or pair-breaking character of 
the individual normal modes. In particular it is found that the half-breathing modes 
give a synergic contribution to the purely electronic pairing; since they are believed to 
be mainly involved in optimally doped cuprates, our findings suggest a joint mechanism, 
with the Hubbard model that captures a crucial part of the physics. 

8. Conclusions and Outlook 

We have presented the following evidence that the W = pairs are the quasi-particles 
that, once dressed, play the role of Cooper pairs: 1) as two-body states they do not feel 
the large on-site repulsion, that would come in first-order perturbation in any theory of 
pairing with any other kind of pairs. 2) The indirect interaction with the background 
particles gives attraction, and bound states with physically appealing binding energies. 
3) The same results are also borne out by exact diagonalization in finite clusters, if 
and only if they have the correct symmetry and filling to give raise to W = pairs in 
partially filled shells. 4) Both in clusters and in the plane, the superconducting flux 
quantization results from the symmetry properties of = pairs. 

The setup of our theory of the effective interaction Wes between two holes (or 
electrons) is quite general; although we developed it in detail for Hubbard Models, it 
is well suited to be extended in order to include other ingredients, hke phonons. In 
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principle we can obtain W^s by our canonical transformation, including systematically 
any kind of virtual intermediate states. The closed-form analytic expression of PVefr 
we obtained includes 4-body virtual states. This describes repeated exchange of an 
electron-hole pair to all orders. One can envisage the pairing mechanism by spin-flip 
exchange diagrams that are enhanced by the 6*4^ symmetry. The argument does not 
depend in any way on perturbation theory, and the equations retain their form, with 
renormalized parameters, at all orders. We find that in the three-band Hubbard Model, 

pairs are more tightly bound close to half filling, but ^-82 pairs are favored when the 
filling increases. Here, like in the BCS theory, the superconducting flux quantization 
relies on the existence of bound pairs of different symmetries. We recall, however, that 
these symmetry labels are not absolute, but depend on the choice of a gauge convention. 
We get attraction and pairing at all fillings we have considered (above half filling), but 
the binding energy of the ^^2 pairs drops by orders of magnitude as the filling increases. 
The three-band Hubbard Model might be too idealized to allow a detailed quantitative 
comparison with experiments; however some qualitative answers are very clear. 

Using the Lieb's theorem |25j and a symmetry analysis based on the W = Q theorem, 
one can fix the way the interactions remove the high degeneracy of the Hubbard Model at 
half filling. This can be used to study pairing in the doped system by the above canonical 
transformation. Exact numerical data on the 4x4 square lattice are available |32j. and 
this allows to test the W = Q pairing mechanism within the one-band Hubbard Model, 
improving over previous weak coupling analvses [48j j49j . For a full application of the 
W=0 theorem, we had to include nontrivial symmetries which are no isometrics; this 
enabled us to classify the W = Q channels and calculate the binding energies analytically. 
Once the analysis is so complete, the criteria for pairing are unambiguous. Moreover 
the success of the weak coupling approach has been noted by other authorsIHZl, and the 
reason is fully clarified by the W = Q theory. 

Despite some evident and important differences that we pointed out above, our 
mechanism can be thought of broadly speaking as a lattice realization of the pioneering 
one proposed long ago^S] by Kohn and Luttinger. The repulsion is avoided by the 
IV = configuration mixing without using parallel spins and high angular momenta. 
The basic source of attraction, however, is essentially the same, since in second-order, in 
the singlet channel, the spin-flip diagram is the only one that survives. We found pairing 
in the singlet channel in a variety of models including carbon nanotubes[72j. This real- 
space approach is suitable for realizing thought experiments, like those involving the 
SFQ (Section El). 

On the other hand, the above results also prove that important ingredients are still 
missing and must be included. The 4x4 model shows evidence of bound pairs of non- 
vanishing momentum, in degenerate representations. This opens up the possibility of 
charge inhomogeneities and Jahn- Teller distorsions. However modeling vibration effects 
in a CUO4 cluster we find that the outcome depends on which phonon mode is most 
strongly coupled. The bound pairs can be flexible enough to prevent distortions and 
actually gain binding energy in the presence of A and B phonons, whilst E vibrations 
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